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Abstract: Shape memory alloys (SMAs) offer a high power-to-weight ratio, large recovery 
strain, and low driving voltages, and have thus attracted considerable research attention. 
The difficulty of controlling SMA actuators arises from their highly nonlinear hysteresis and 
temperature dependence. This paper describes a combination of self-sensing and 
model-based control, where the model includes both the major and minor hysteresis loops as 
well as the thermodynamics effects. The self-sensing algorithm uses only the power width 
modulation (PWM) signal and requires no heavy equipment. The method can achieve 
high-accuracy servo control and is especially suitable for miniaturized applications. 
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1. Introduction 

Shape memory alloys (SMAs) are metals that can recover from strains of up to 10% via 
stress- or temperature-induced crystalline transformation between high-temperature austenite and 
low-temperature martensite phases [1]. The SMA actuator is constructed from a fiber-like SMA wire 
designed to contract and extend like real muscles. At room temperature, an SMA wire is soft and pliable, 
very much like a nylon thread. However, when heated it begins to contract sharply with a large force 
and eventually becomes as stiff as a piano wire. The maximum strain is typically 4.5% of its original 
length. When the SMA is again cooled to room temperature it softens and recovers its original length. 
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Due to the characteristics of a high power-to- weight ratio, large recovery strain, and low driving 
voltages, the SMA actuator has been used in wide variety of applications including aircraft wing 
controls [2,3], robotic grippers [4-7], automotive mirror actuators [8], active vibration suppression [9], 
active endoscopes [10], and legged robots [11,12]. Most SMA applications require some form of length 
control, and a simple implementation involves using separate strain sensors for the SMA deformation 
for feedback control; however, this can be very difficult for some miniature applications and the 
sensorless approach offers an attractive alternative. The sensorless SMA control appraoches can be 
divided into two major categories. The first approach is to use the so-called self-sensing properties of 
the SMA actuator, whereby the change in the SMA electric resistance is measured to estimate the 
corresponding strain. Curve-fitting and a neural network have been used to model the SMA self- sensing 
properties [5,13]. These models were able to describe the major hysteresis loop of the SMA actuator 
but not the minor hysteresis loops. Most of the control applications also employed conventional PD 
control for the feedback action. The second approach to the sensorless SMA control uses no 
measurement feedback, but depends instead on mathematical models to estimate the SMA 
strain [14-16]; obviously this method is sensitive to the accuracy of the mathematical models. 

In this paper, we propose a modified approach for precision sensorless SMA servo control that 
consists of three components: (1) a hysteresis model that combines the strengths of the two sensorless 
control strategies, (2) a thermodynamics model to compensate for the temperature effect, and (3) a 
spring model to include the strain energy effect. The hysteresis model is based on the Duhem differential 
model, and is used to describe both the major and minor hysteresis loops. A detailed model is necessary 
to fulfill the stringent precision control requirements. Variable supply voltages have previously been 
used to induce the SMA self-sensing relationship [5,13]. However, the resulting device (i.e., the 
programmable power supply) is far too large for miniature robot applications; this difficulty is overcome 
in this paper by using the PWM signal to obtain the self-sensing characteristics. This requires only a 
constant supply voltage, and so a bulky programmable power supply can be replaced by very small IC 
chips that are more suited for miniature applications. The PWM signal is also employed for controlling 
the SMA actuator displacements. 

The Section 2 of the paper describes the experimental setup for the proposed control scheme. 
Section 3 explains the self- sensing property of the SMA actuator. Modeling of the SMA actuator is 
described in Section 4. The proposed scheme of tracking control based on self-sensing feedback and 
inverse hysteresis compensator is described in Section 5. Section 6 shows the experimental results. 
Finally, the conclusions are given in Section 7. 

2. Experimental Setup 

The main components of the experimental setup (a test platform and an electric circuit) are shown in 
Figure 1. The test platform was used to investigate the characteristics of a bias-type SMA wire actuator 
(the SMA wire contracts when heated, and it expands with the aid of the bias spring when cooled). 
A 167-mm-long NiTi-based SMA actuator with a diameter of 150 jam (BMF150, TOKI) was installed 
on the test platform. One of the ends of the SMA actuator was fixed to the platform while the other end 
was connected to a moving slider. A pair of linear guides restricted the slider to move only in one 



Sensors 2010, 10 



114 



dimension horizontally. The linear bias spring provided a restoring force to the SMA actuator. In 
addition, an LED displacement sensor (OMRON) was integrated to measure the displacement of the 
SMA actuator. Note that the displacement sensor was used in this study simply to validate the control 
results, and not to provide a feedback signal to the controller. 

Figure 1. Experimental setup. 




Linear guide 

A schematic of the electric circuit is shown in Figure 2. A multifunction data acquisition card 
(+5 V full-scale range, 12-bit resolution; PCI-1711, Advantech) was employed to send the PWM signal 
via the digital output and measure the amplified voltage V amp via the analog input. A Darlington driver 
IC (ULN2003AP) was used as a switching element to control the heating or cooling state of the SMA 
actuator. A DC voltage source V s was connected to the SMA actuator to supply a DC voltage to heat 
the SMA actuator. A resistor, R, was connected serially to the SMA actuator to prevent it from 
overheating. VCE of the Darlington driver was then amplified by a differential amplifier to enlarge its 
variation due to the electric-resistance variation of the SMA actuator during the phase transformation 
process. The amplified voltage Vamp was measured by the data acquisition card. 



Figure 2. Schematic of the electric circuit. 
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3. Self-Sensing Property of SMA Actuator 



Figure 3 shows the relationship between the PWM signal that inputs to the Darlington driver and the 
measured V amp . The supply voltage V s is set to 6 V and 100 Hz PWM signal with 40% duty ratio is input 
to the Darlington driver. It can clearly be seen in Figure 3 that during the "turn-off' state (low voltage 
level) of the input PWM signal, the corresponding V amp is saturated at 5 V. On the other hand, during 
the "turn-on" state (high voltage level) of the input PWM signal, the corresponding V amp drops to some 
steady value (in this case about 0.7 V) after 1ms. At each falling edge of the input PWM signal, a V amp 
value is recorded as Va- In other words, one Va is recorded during each PWM cycle. Note that it takes 
at least 1ms for the V ump to drop from 5 V to a steady value; therefore, the duty ratio of the 
input 100 Hz PWM signal must higher than 10% in order to measurement the correct Va- 

Figure 3 shows the relationship between the PWM signal supplied to the Darlington driver and the 
measured V amp . The supply voltage V s was set to 6 V and a 100 Hz PWM signal with a 40% duty cycle 
was supplied to the Darlington driver. The figure clearly shows that V amp saturated at 5 V during the 
"turn-off state (low voltage level) of the input PWM signal, and dropped to some steady value (in this 
case about 0.7 V) after 1 ms during the "turn-on" state (high voltage level). The V amp value at each 
falling edge of the input PWM signal is recorded as Va', that is, one Va value is recorded during each 
PWM cycle. Note that it took at least 1 ms for the Vamp to drop from 5 V to a steady value; therefore, 
the duty cycle of the input 100-Hz PWM signal had to be higher than 10% to allow the correct Va to 
be measured. 

Figure 3. Relationship between V amp (blue solid line), input PWM (red dashed line) and Va (circle). 
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An open-loop experiment was performed to investigate the relationship between the contraction 
displacement of the SMA actuator (D s ) and V A - The room temperature was 23 °C and no additional 
cooling method was applied to the test platform. At the beginning of the experiment, the 100 Hz PWM 
signal with a 40% duty cycle was used as the input signal to slowly heat the SMA actuator. After the 
SMA actuator was fully contracted, the duty cycle of the input PWM signal was changed to 15%. The 
SMA actuator was then cooled gradually. The values of D s and Va were recorded throughout the 
heating and cooling processes. The experiment was repeated 10 times, and Figure 4 shows the resulting 
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plot of D s versus V A . The data show high repeatability and are modeled in Section VI so as to 
characterize the self-sensing relationship of the SMA actuator. Note that the D. s - versus- Va plot varies 
with the pretension of the bias spring [13]. In this study, the pretension was set to 65 gw, and the spring 
constant (k s ) of the bias spring was 35.035 N/m. 

Figure 4. Plot of D s versus V A for 10 experiments. 




4. Modeling of SMA Actuator 



4.1. Modeling of Hysteresis 



An SMA actuator exhibits the hysteresis phenomenon [1]. Figrue 5 shows the relationship between 
temperature and strain of the BMF150 device as published in the product specifications [17]. For the 
heating/contraction process, the strain was small until the temperature was above 75 °C, whereas for the 
cooling/elongation process a large strain appeared after the temperature decreased to below 70 °C. 

Figure 5. Relationship between strain and temperature for BMF150 [17]. 
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The Duhem differential model [16,18] is used to model the hysteresis phenomenon of the 
SMA actuator: 



(1) 



I y(t) = g + (u(t), y(t))(u(t)) + ~ g-(u(t), y(t))(u(t)T 

i y(0) = y 0 , 

where u(t) denotes the input, y(t) is the output, y 0 is the initial value of the output, g ± is slope function, 
the subscripts + and - represent increasing and decreasing curves, respectively, and: 

(u(t)f=(\u(t)\±u(t))/2 (2) 

Equation (2) indicates that (u(t)) + = u(t) and (u(t)) =0 if u(t)>0 , and (u(t)) + =0 and 
(ii(0)~ = u(t) if u(t) > 0 . As a result, the slope function of (1) is g+ when u(t) > 0 and g. when u(t) < 0 . 

The Gaussian combination membership function (Gcmf) is chosen as the slope function of the 
hysteresis model. The G C mf is a combination of two Gaussian membership functions: 

G CMf( U ) = k CMF ■ G MF,M) ■ G MF,2^ ( 3 ) 

where k C MF is the gain, G MF l (u) and G MF 2 (u) are the modified Gaussian membership functions, and is 
defined as: 



G MF,\ (") : 



MF,2 



(«) = 



(l-c x )exp 
1, 
1, 

(l-c 2 )exp 



2of 



(u-ju 2 ) 



+ q , if u < /u { 
ifu> J u 1 

if u < ju 2 
+ c 2 , if u > ju 2 



(4) 



(5) 



J 



where //, denotes the mean, o i denotes the variance, and c, is an offset value. Figure 6 gives an example 
of a Gcmf based on two Gaussian membership functions: G MF [With fx\ = 20, <jf= 9, c\ = 0.1, and G MF2 
with // 2 = 60, CTj = 100, c 2 = 0.1. The hysteresis between strain s and temperature T shown in Figure 5 is 
modeled by the differential equation: 



d^\G CMF ^T), ifj>0 
dT \G CMF _(n ifT<0 



(6) 



where G CMF + and G CMF _ are the slope functions for heating and cooling curves, respectively. Figure 7 

compares the hysteresis phenomenon of the SMA actuator between the device specifications and the 
model given by (6). The corresponding modeling parameters are listed in Table 1. 
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Figure 6. Example of a Gaussian combination membership function (dotted line) based on 
two modified Gaussian membership functions G M fj (solid line) and G M f,i (dashed line). 
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Figure 7. Comparison of the hysteresis of the SMA actuator between the device 
specifications (solid line) and the modeling result (dashed line). 
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Table 1. Parameters for SMA hysteresis modeling. 
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To model the minor hysteresis loops, (6) is modified by multiplying the slope function by a gain: 

h_(T)-£ „ 



ds 
dT 



h(T)-h + (T) 
s-K(T) 



'CMF,+ 



(T), if T > 0 



(7) 



'CMF, 



XT), ifr<o 



h_(T)-h + (T) 

where h + (T) and h_(T) are the increasing and decreasing curves of the major hysteresis loop computed 

by (6), respectively The simulated hysteresis model of the SMA actuator based on (7) is shown in 
Figure 8. 

Figure 8. Simulated hysteresis characteristics of SMA actuator. 
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4.2. Modeling of Inverse Hysteresis 



An inverse hysteresis model has to be derived to compensate for the hysteresis of the SMA 
actuator. We derived the following inverse hysteresis model by inverting the hysteresis differential 
equation (Equation 7) [19]: 

'h_(T)-h(T) 1 



dT_ 

ds 



h_(T)s + S G CMF+ (T) + S 
h_(T)-h + (T) 1 
e-h + {T) + S G CMF ^T) + S' 



if s>0 
if s<0 



(8) 



where d is a positive arbitrarily small constant that allows the differential equation to be solved. Figure 9 
illustrates the simulated inverse hysteresis model of the SMA actuator based on (8). 
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Figure 9. Simulated inverse hysteresis characteristics of SMA actuator. 




4.3. Modeling of Temperature Dynamics 



The temperature dynamics are modeled by the heat transfer equation, which balances the heat across 
the SMA actuator. According to the Joule effect, the SMA actuator is heated when electric current is 
passed through it. We assume that heat loss occurs only via natural convection. The temperature 
dynamics are given by the following differential equation [14,16,20]: 

T = —{P T -xd s L 0 h{T-T amb )) (9) 
mc 

where T amb is the ambient temperature; T, m, c, d s , and L 0 are the temperature, mass, specific heat, 
diameter, and undeformed length of the SMA actuator, respectively; P T is the electric power required to 
change the temperature of the SMA actuator; and h is the coefficient of convectional heat transfer. Note 
that h for a thin metal wire under natural cooling varies with temperature [21], and so its 
value was obtained by the method provided in [21] under the following conditions: T am b = 23 °C, 
d s = 1.5 x 10"* m, L 0 = 0.167 m, and an air pressure of 1 atm. Using the MATLAB curve-fitting tool, 
the relationship between heat convection coefficient h and temperature T is represented as: 

h(T) = ai e aiT +b x e hT (10) 

where a x = 85.28, a 2 = 0.001727, b x = -106.4, and b 2 = -0.08706. In addition, the specific heat of the 
SMA actuator does not remain constant during the actuation process, instead differing between the 
martensite phase (c m ) and the austenite phase (c a ); accordingly, the value of c is based on the percentage 
of martensite transformation: 

c = j-c M +(L~)c t (11) 



where L 0 is the undeformed length of the SMA actuator and L is the actual length of the SMA actuator 
(L = L 0 -D s ). 
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4.4. Electric Power Calculation 

The electric power provided by the power source not only heats the bias-type SMA actuator but also 
supplies energy to elongate the bias spring. The total power, P, is given by: 

P = P T +P W (12) 

where P T is the electric power required to change the temperature of the SMA actuator, as given in (9), 
and P w is the work required to provide the elastic potential energy to the bias spring, which is derived 
from the potential energy of the bias spring according to: 

W = ±D s [k s l 0 + k s (l 0+ D 5 )] 

where k s is the spring constant of the bias-spring, D s is the displacement of the SMA actuator, and l 0 is 
the initial length of the bias-spring. P w is given by: 

P w =W = k s (D s +l 0 )D s (13) 
The duty cycle of the PWM signal D that is supplied to the SMA actuator can then be estimate as 

D 'lv^r- (14) 

V PWM PWM 

where V PW m is the high-level voltage of the PWM signal and I PWM is the electric current across the SMA 
actuator when D = 100%. 

4.5. Modeling of Self- Sensing Properties 

To model the self-sensing properties of the SMA actuator, the curve of D s versus V A shown in 
Figure 4 is separated into three phases. In phases I and II, the mathematical models are obtained using 
the curve-fitting technique. The two heating curves in phases I and II, and the cooling curve in phase II 



are modeled respectively by the following equations: 

D s = f(V A ) = -0.8554V, 2 + 2.7850V, -0.9077 (15) 

D s = f(V A ) = 0.0192V] +0.0081V A 2 +1.2960V, -0.3065 (16) 

D s = /(V A ) = 0.0811^+1.3160^ -0.4681 (17) 



The D s - versus- V A curve in Figure 4 exhibits slight hysteresis. However, the mathematical models 
given in (15)— (17) only describe the major hysteresis loop. The Duhem differential model is therefore 
employed to obtain the complete model including both major and minor hysteresis loops. The 
self- sensing properties of the SMA actuator are given by: 
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dV A 



h_(V A )-h + (V A ) 

d s -K(v a ) 



S + (V A ), ifV A >0 
S_(V A ), HV A <0 



(18) 



[h_(v A )-h + (y A ) 

where S + (V A ) and S_(V A ) are the slope functions of the increasing and decreasing curves, respectively, 
which can be obtained by differentiating (15)— (17) once. h + (V A ) and h_(V A ) are the increasing and 

decreasing curves, respectively, of the major loop computed by the following equation: 

dp, = jXcrj, ify A >o 
dV A \s_(v A ). av A <o 

Figure 10 shows the mathematical model comparing with the experimental data. 



(19) 



Figure 10. Comparison of simulation and experimental results for the D s - versus- Va curve. 
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5. Tracking Control with Self-Sensing Feedback and Inverse Compensation 



The SMA actuator is a highly nonlinear system owing to its hysteresis characteristics; accordingly, a 
feedforward inverse compensator is designed to compensate the hysteresis. A schematic of the inverse 
compensator for controlling the displacement of the SMA actuator is depicted in Figure 1 1 . The input 
of the inverse compensator is the reference strain, e, where s = D s I L 0 . Block I outputs the 
corresponding temperature of the SMA actuator to block II, which computes the power (P T ) required 
to heat (or cool) the SMA actuator to the specific temperature. Block III computes the power (P w ) 
required to deform the bias spring. Total power P (= P T + P w ) is supplied to block IV to calculate the 
corresponding duty cycle D of the input PWM signal. In summary, the feedforward inverse 
compensator estimates the duty cycle of the input PWM signal required to heat (or cool) the SMA 
actuator to the desired length. 
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Figure 11. Schematic block diagram of the feed-forward inverse compensator for SMA 
length control. 
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Figure 12 presents the control flow diagram of the proposed self-sensing feedback control with 
inverse hysteresis compensation. The strain estimated by the self-sensing characteristics described in (18) 
is used as the feedback source to compare with the reference strain. The error is supplied to a 
conventional PID controller to generate the appropriate duty cycle, D 2 , of the PWM signal to implement 
the tracking ability. In addition, the reference strain is supplied to the feedforward inverse compensator 
to estimate the duty cycle, D u required for the reference strain. The SMA actuator is controlled by the 
PWM signal with duty cycle D computed by adding D\ and D 2 . 



Figure 12. Flow diagram of the SMA actuator length control. 
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6. Experimental Results 



Two experiments were performed to examine the tracking performance of the proposed control 
architecture. The first experiment used a sinusoidal reference signal. Figure 13 presents the experimental 
measurements from the displacement sensor, showing the reference sinusoidal signal and the results 
obtained using only the PID controller and simultaneously using both the PID controller and the 
feedforward inverse hysteresis compensator. The parameters of the PID controller were identical in both 
cases. Figure 13 shows clearly that the tracking error was large for the PID controller only when the 
reference signal switches between heating and cooling. This was attributed to by the very different 
heating and cooling hysteresis curves. Such errors reduced when the feedforward inverse hysteresis 
compensator was incorporated with the PID controller since this compensates the hysteresis of the 
SMA actuator. The RMS values of the tracking errors are listed in Table 2. 
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The second experiment used a multistep reference signal. The experimental results measured by the 
displacement sensor are shown in Figure 14, which indicates that when using only the PID controller the 
displacement trajectory could not follow the reference well at the beginning (0-2 mm) of the experiment. 
This was due to the nonlinearity of the SMA actuator shown in Figure 5 the temperature of the SMA 
actuator needed to increase about 65 °C to move from 0 to 1% strain, but it only needed an increase 
of 35 °C for the SMA actuator to move from 1% to about 5.7% strain. However, the tracking errors 
improved when the feedforward inverse hysteresis compensator was incorporated with the PID 
controller. Note that the parameters of the PID controller were identical when using the PID controller 
alone or the combination of the PID controller and the feedforward inverse hysteresis compensator. The 
RMS values of the tracking errors are listed in Table 2. 

Figure 14. Experimental results for a multi-step reference input. 
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As a comparison, Figure 15 showed the same control with different models: the Duhem model, the 
major loop model proposed in [13], and the single curve model proposed in [22]. The major loop model 
merely modeled the major hysteresis loop and thus resulted in rapid chattering in the control effort when 
the actuator response approaches steady state. This was due to the fact that the underlying model 
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switched between the "heating" and the "cooling" curves. To avoid this difficulty, the single loop model 
used only a single polynomial for the self-sensing characteristics. The Duhem differential model also 
suppressed the chattering behavior. The RMS tracking error for the major loop model, the single curve 
model, and the Duhem model were 0.2890, 0.2781, and 0.2311 mm, as shown in Table 2. It could be 
seen that the control based on the Duhem model and teh single curve model exhibited larger transient 
responses than the control with the proposed model, which achieved an RMS error of 0.1223 mm. 

Figure 15. PID control with different self-sensing models: the single curve model 
(dash-dotted line), the major loop model (dashed line), and the Duhem model (solid line). 
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Table 2. RMS errors of the tracking errors. 
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7. Conclusions and Future Remarks 

This paper described a PWM based self- sensing feedback controller with inverse hysteresis 
compensator for a SMA actuator. The proposed SMA compensator comprised an inverse hysteresis 
model to represent the major and minor hysteresis loops, a temperature dynamics model to compute the 
required input power to heat up the SMA actuator, and a spring force model that took accounted for the 
strain energy required to deform the actuator. The inverse hysteresis model was based on the Duhem 
self-sensing characteristic represented by the D S -V A relationship. Both the major hysteresis loop and the 
minor hysteresis loops were considered. As a result, the model enabled accurate estimate of the actuator 
strain by using the electrical potential across the actuator. Experimental results showed that the 
self-sensing model achieved smaller transient error and can effectively suppress the 
chattering phenomenon. 

It is worth noting that the current control scheme depends on the precise knowledge of the ambient 
temperature and the material properties. While the material properties remain mostly constants, 
the ambient temperature may experience unexpected changes. Future research will address the 
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robustness issues of the proposed control against unknown ambient temperature changes and 
material property variations. 
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